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I.  INTRODUCTION 


The  use  of  hydrogen  masers  as  frequency  standards  on  spacecraft  requires 
careful  consideration  of  long-term  reliability  issues.  A  primary  issue  is  the 
maser's  hydrogen  budget,  because  many  of  the  common  failure  modes  of  a 
hydrogen  maser  involve  either  the  hydrogen  source  (storage  or  dissociator)  or 
sink  (ion  pump  or  getters).  An  excess  of  hydrogen  atoms  in  the  F  =  1,  M  =  0 
state  must  be  continuously  fed  into  the  maser  bulb.  The  ratio  of  the  rates  at 
which  atoms  in  the  proper  hyperfine  state  enter  the  bulb  and  molecules  are  fed 
into  the  dissociator  will  impact  strongly  the  maser's  long-term  reliability. 
This  ratio  will  be  determined  mainly  by  the  dissociator' s  efficiency  and  the 
properties  of  the  state-selecting  magnet.  Atomic  state  selection  is  performed 
in  hydrogen  masers  by  a  quadrupole  or  hexapole  magnet,  which  focuses  atoms  in 
the  F  =  1,  M  =  0  state  at  the  bulb's  entrance  orifice  and  defocuses  the  F  =  0, 
M  =  0  atoms.  The  focal  length  of  those  magnets  depends  strongly  on  atomic 
speed,  and  typically  only  a  narrow  range  of  atomic  velocities  will  be 
focused.  Clearly,  dissociators  yielding  narrower  velocity  distributions  could 
more  efficiently  use  the  hydrogen  supply. 

Relatively  little  is  known  of  the  velocity  distribution  of  atoms  effusing 
out  of  radio  frequency  (rf)  discharge  hydrogen  dissociators.  In  many  cases, 
the  tacit  assumption  is  made  that  the  atoms  will  be  in  thermal  equilibrium 
with  the  dissociator  wall,  but  that  is  not  necessarily  true.  The  threshold 
for  molecular  dissociation  by  collisions  with  discharge  electrons  is  about  8.5 
eV.  Because  the  molecular  binding  energy  is  approximately  4.7  eV,  each  atom 
carries  away  an  excess  energy  of  about  2  eV  (Ref.  1).  Depending  on 
dissociator  geometry  and  gas  density,  the  hydrogen  atoms  may  or  may  not 
undergo  enough  gas  and  wall  collisions  to  thermal ize  fully.  Walraven  and 
Silvera  (Ref.  2)  have  studied  the  characteristics  of  a  beam  of  hydrogen  atoms 
produced  in  a  microwave  discharge  followed  by  a  thermal  accommodator  and  found 
the  velocity  distributions  to  be  Maxwellian  at  the  accommodator  temperature. 
Hershcovitch  et  al .  (Ref.  3)  studied  a  similar  dissociator-and-accommodator 
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combination,  operated  at  higher  pressures  and  flow  rates,  and  measured 
velocity  distributions  in  reasonable  agreement  with  the  calculated  supersonic- 
flow  distributions.  Because  thermal  accommodators  are  not  commonly  used  in 
maser  dissociators ,  these  results  are  not  directly  applicable  to  our 
problem.  Miller  (Ref.  4)  investigated  a  helium-cooled  microwave  discharge 
dissociator  at  relatively  high  pressure  (3.3  Torr)  and  measured  a  velocity 
distribution  that  was  slightly  narrower  than  the  Maxwellian  at  the  estimated 
discharge  temperature.  We  are  studying  the  velocity  distributions  of  atomic 
hydrogen  produced  in  an  rf  discharge  dissociator  having  a  geometry  and 
operating  parameters  resembling  those  of  a  maser  dissociator.  The 
experimental  technique  and  preliminary  results  are  presented  in  this  report. 
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II.  EXPERIMENTAL  APPARATUS 


Atomic  hydrogen  velocity  distributions  have  been  determined  using  a 
magnetic  deflection  technique.  Figure  1  exhibits  the  experimental 
arrangement.  Hydrogen  is  fed  through  a  standard,  temperature-controlled  Pd-Ag 
alloy  leak  (Ref.  5)  into  a  cylindrical,  double-walled  Pyrex  bulb.  The  bulb  is 
15  cm  long,  with  a  1.9  cm  internal  diameter  and  a  3.8  cm  external  diameter. 
Compressed  air  flows  between  the  walls  to  provide  cooling,  and  rf  power  is 
inductively  coupled  to  the  discharge  through  an  external  25-turn  coil.  The 
hydrogen  beam  exits  the  dissociator  through  a  slit  that  is  0. 1  cm  long  and 
0.025  cm  wide,  and  is  collimated  by  a  second  slit  that  is  0.025  cm  wide,  set 
at  d  =  63.7  cm  away  from  the  source  slit.  The  beam  then  travels  between  the 
pole  pieces  of  an  L  =  11.4  cm  long  electromagnet  configured  in  the  "two-wire" 
geometry  described  by  Rabi  et  al.  (Ref.  6).  After  additional  travel  through  a 
drift  tube  that  is  D  =  71.3  cm  long,  the  beam  is  detected  by  a  quadrupole  mass 
analyzer . 

Molecular  hydrogen,  having  no  magnetic  dipole  moment,  will  travel  through 
the  magnet  without  deflection,  but  the  hydrogen  atoms  will  be  deflected  by  the 
inhomogeneous  magnetic  field.  The  dissociator  is  attached  to  the  rest  of  the 
apparatus  by  flexible  vacuum  bellows  and  can  be  displaced  transversally  by 
micrometer  screws.  In  this  way,  the  angular  distribution  of  atoms  deflected 
by  the  field  can  be  measured. 
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HYDROGEN  ‘  TWO-WIRE”  QUADRUPOLE  MASS 

DISSOCIATOR  MAGNET  FILTER  AND  DETECTOR 


Fig.  1.  Schematic  View  of  the  Experimental  Arrangement. 

The  hydrogen  dissociator  can  be  displaced  transversal ly . 
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III.  ATOMIC  DEFLECTION  IN  A  "TWO-WIRE"  MAGNETIC  FIELD 


An  atom  having  magnetic  moment  u,  immersed  in  a  magnetic  field  fi,  has  an 
energy  E  =  The  force  acting  on  the  atom  is  f  =  -grad(y-fi)  *  grad 
(H),  where  is  the  component  of  the  atomic  magnetic  moment  along  the  field 
direction.  In  general,  will  depend  on  the  atomic  ground  state  quantum 
numbers  and  the  magnetic  field  strength.  If  an  atom  has  total  electron 
angular  momentum  J  =  1/2  and  nuclear  spin  I,  both  in  units  of  h/2it,  then  its 
effective  magnetic  moment  is  given  by  Breit-Rabi's  formula  (Ref.  7) 


peff  =  1 


e  +  2M/(2I  +  1) _ 

[1  +  4Me/(2I  +  1)  «•  e2]1/2 


(1) 


where  the  sign  corresponds  to  total  angular  momentum  F  =  I  -  1/2,  the 
sign  to  F  =  I  +  1/2,  and  M  is  the  total  azimuthal  quantum  number.  e  is 
proportional  to  the  ratio  of  magnetic-to-hyperfine  energy,  e  =  gvQH/W,  where  g 
is  the  Lande  factor  for  the  atom,  Wq  is  the  Bohr  magneton,  and  W  is  the  atomic 
hyperfine  splitting. 

Figure  2  exhibits  the  cross  section  of  the  "two-wire"  magnet;  2a  is  the 
separation  between  the  "equivalent  wires"  (Ref.  8).  If  a  narrow  atomic  beam 
enters  the  magnet  at  x  =  1.2a  (as  determined  by  the  collimating  slit)  and  the 
atomic  transversal  displacements  within  the  magnet  are  small,  then  the  force 
acting  on  the  atoms  remains  approximately  constant  (Ref.  9),  and  so  does  the 
acceleration : 


a 


K 


0.984  weff 
a  m 


H  i 


(2) 


where  H  is  t|ie  magnetic  field  intensity  at  the  atomic  beam,  m  is  the  atomic 
mass,  and  i  is  a  unit  vector  pointing  along  the  x  axis.  Under  these 
conditions,  the  trajectory  equations  for  the  atoms  traveling  from  source  to 
detector  can  be  solved  easily.  With  no  magnetic  field,  an  atom  of  magnetic 
moment  and  speed  V  leaving  the  source  slit  at  xs  will  be  detected  at  xD 
=  -Xg  (L  -*•  D)/d.  If  the  magnetic  field  is  set  at  H  and  the  detector  is  not 


9 


BEAM 


moved,  an  identical  atom  will  have  to  leave  the  plane  of  the  source  slit  at 
xs' ,  given  by 


x 


s 


_D_1 

L  +  DJ 


(3) 


in  order  to  be  detected.  If  fQ(xg)  is  the  distribution  of  detected  atoms  as  a 
function  of  source  slit  position  at  zero  magnetic  field,  and  g ( V )  is  the  speed 
distribution  of  atoms  leaving  the  source  slit,  then  the  distribution  of 
detected  atoms  at  nonzero  field,  f(xs),  will  be  given  by 

a> 

f(xs)  *  I  f  f0  (xs  -  Kfm/V2)  g( V)  dV  (4) 

F ,  M  o 

where 


.  0-492 
FM  a 


yeff 

m 


H  Ld  (1 


+ 


i-] 

L  +  DJ 


(5) 


Both  fQ(xs)  and  03,1  be  ®easured»  and  then  eQ*  (*0  can  be  used  to  verify 

whether  a  given  velocity  distribution  is  consistent  with  these  measured 
distributions. 
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IV.  CALIBRATION  AND  PROOF-OF -CONCEPT  EXPERIMENT 


A  preliminary  experiment  to  study  the  deflection  of  a  rubidium  beam  was 
performed  to  calibrate  the  electromagnet  and  test  the  use  of  Eq.  (M)  to  verify 
an  assumed  velocity  distribution. 

Equation  (1)  reveals  that  for  an  atom  having  nuclear  spin  I,  =  0  if 
e  =  -2M/(2I  +1).  At  field  intensities  determined  by 


-M 

=  21  +  1 


W_ 

v0 


(6) 


atoms  with  M  <  0  will  not  be  deflected  by  the  field.  Natural  rubidium  is  a 
mixture  of  72.2 %  Rbgg  (I  =  5/2)  and  27.8)1  Rbgy  (I  =  3/2).  Rbgg  will  have  zero 
effective  moments  at  H  =  361  G  (M  =  -1)  and  H  =  722  G  (M  =  -2).  Rbg^  will 
have  a  zero  effective  moment  at  H  =  1221  G  (M  =  -1).  Figure  3  exhibits  the 
detected  atomir  .signal  as  a  function  of  magnet  coil  current,  with  the  rubidium 
beam  source  and  detector  on  the  apparatus  axis.  At  the  field  values  just 
listed,  atoms  in  the  zero  effective  moment  states  will  travel  through  the 
magnet  without  deflection,  leading  to  an  increased  signal.  The  effect  of  the 
zero  effective  moment  states  can  be  made  more  apparent  on  the  figure  by 
subtracting  the  local  sloping  background  from  each  one  of  the  marked 
features.  The  peaks  labeled  A  and  B  can  then  be  identified  with  the  zero 
effective  moments  at  722  and  1221  G,  respectively,  yielding  a  calibration  of 
66.3  G/A  for  the  magnet. 

The  rubidium  source  for  these  measurements  was  a  two-chambered  oven 
designed  and  operated  to  ensure  effusive  flow.  Under  those  conditions,  the 
atomic  beam  speed  distribution  should  be  beam-Maxwellian: 


g(V)  =  (2/e)  (V/8)3  exp  [-(V/B)2] 


(7) 


where  6  =  /(2kT/m)  and  T  is  the  oven  temperature.  We  measured  f(xs)  for  H  = 
1061  G,  H  =  1220  G,  and  H  =  1525  G,  and  also  measured  fQ  (xg).  We  then  used 
Eq.  (4),  with  the  speed  distribution  given  by  Eq.  (7),  to  calculate  the 
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Fig.  3 


COIL  CURRENT  (A) 

Rubidium  Atomic  Beam  Signal  (Detector  and  Oven  on  Axis) 
vs  Magnet  Coil  Current.  Arrows  exhibit  features  resulting 
from  zero  effective  moment  states.  Dashed  lines  were 
obtained  by  subtracting  the  local  constant-slope 
background  from  each  feature.  A:  Rbgc,  M  =  -2; 

B:  Rbg7,  M  *  -1. 
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expected  f(Xg)  at  those  field  values.  The  results  are  presented  in  Fig.  4. 
For  x  <  0,  the  agreement  between  calculated  and  measured  detected  atom 
distributions  is  somewhat  poor.  This  side  of  the  distribution  is  near  the 
convex  pole  piece  of  the  magnet,  and  the  constant-force  approximation  made  to 
derive  Eq.  (3)  breaks  down  for  atoms  traveling  too  close  to  the  pole  piece. 
For  x  >  0,  the  agreement  between  calculated  and  measured  f(x£)  is  excellent. 
These  results  fully  validate  our  experimental  technique. 
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SIGNAL 
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0  -0.2  0  0.2  -0.2  0  0.2  0.3 

DETECTOR  POSITION  (cm) 


Fig.  4.  Detected  Rubidium  Atom  Distributions  at  Indicated  Magnetic 
Fields.  Oven  temperature:  484  K.  Solid  line:  measured 
distributions.  Dashed  line:  calculated  distributions  for 
a  Maxwellian  beam. 
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V.  ATOMIC  HYDROGEN  VELOCITY  DISTRIBUTIONS 


Figure  5  reveals  the  atomic  hydrogen  beam  fraction  as  a  function  of  total 
pressure  in  the  dissociator  bulb,  and  discharge  rf  power.  These  data  were 
obtained  at  a  discharge  frequency  of  205  MHz.  The  discharge  power  was 
corrected  for  reflection,  and  the  atomic  fraction,  measured  with  the 
quadrupole  mass  analyzer,  was  corrected  for  the  H+  background  when  detecting 

h2. 

We  have  taken  preliminary  atomic  hydrogen  deflection  data  at  total 
dissociator  pressures  of  0.2,  0.4,  and  0.7  Torr,  and  at  995  and  1326  G 
magnetic  fields.  Zero-field-detected  atom  distributions  have  also  been 
measured  at  each  operating  pressure.  We  then  used  Eq.  (4)  with  the  bearo- 
Maxwellian  distribution  given  by  Eq.  (7)  to  test  the  validity  of  the 
hypothesis  of  a  thermal  beam.  The  source  temperature  T  was  treated  as  a  free 
parameter.  Figure  6  exhibits,  as  an  example,  the  deflection  data  at  995  G  and 
0.4  Torr,  together  with  the  expected  detected  atom  distributions  for 
Maxwellian  beams  at  3 00  and  5 00  K.  As  shown  by  Eq.  (3),  the  transversal 
atomic  displacement  by  the  field  is  inversely  proportional  to  V2.  Therefore, 
it  is  clear  that  the  distribution  at  300  K  contains  a  large  excess  of  slow 
atoms  (large  deflections),  while  the  distribution  at  500  K  contains  a 
correspondingly  large  excess  of  fast  atoms  (small  deflections).  Similar 
problems  were  encountered  when  analyzing  the  rest  of  the  deflection  data. 
These  results  indicate  that  the  atomic  velocity  distribution  in  the  beam  must 
be  significantly  narrower  than  the  beam-Maxwellian . 

To  estimate  qualitatively  the  characteristics  of  the  actual  velocity 
distribution  in  our  atomic  hydrogen  beam,  we  decided  to  model  it  by  a  Gaussian 
distribution  of  peak  speed  Vq  and  width  parameter  o: 

g(V)  =  ^  exp  (-(V  -  V0)2/2o2]  (8) 

Vq  and  o  were  treated  as  free  parameters  and  adjusted  to  provide  a  reasonable 
approximation  to  the  measured  deflection  data.  Given  the  preliminary  nature 
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DISSOCIATION  FRACTION 


Fig.  5.  Atomic  Hydrogen  Beam  Fraction  vs  rf  Discharge 

Power.  Total  dissociator  bulb  pressures  are  indicated. 
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NORMALIZED  SIGNAL 


SOURCE  POSITION  (cm) 


Fig.  6.  Detected  Hydrogen  Atom  Distribution  at  995  G  and  0.40  Torr. 
Dots:  measured  distribution.  Curves:  calculated 
distributions. 
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of  the  deflection  data  and  model  velocity  distribution,  no  attempt  was  made  to 
obtain  a  "best  fit."  Figures  7,  8,  and  9  exhibit  the  deflection  data  at  995  G 
for  dissociator  total  pressures  of  0.2,  0.4,  and  0.7  Torr,  with  the 
corresponding  results  calculated  using  the  velocity  distributions  given  by  Eq. 
(8).  A  much  closer  agreement  with  the  measured  data  is  obtained  than  when 
using  Maxwellians.  Similar  results  were  obtained  for  the  data  at  1326  G. 
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NORMALIZED  SIGNAL 


Fig.  7.  Detected  Hydrogen  Atom  Distribution  at  995  G  and  0.20 
Torr.  Dots:  measured  distribution.  Curve: 
calculated  distribution. 
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NORMALIZED  SIGNAL 


Fig.  8.  Detected  Hydrogen  Atom  Distribution  at  995  G  and  0.40 
Torr.  Dots:  measured  distribution.  Curve: 
calculated  distribution. 
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NORMALIZED  SIGNAL 


Fig.  9.  Detected  Hydrogen  Atom  Distribution  at  995  G  and  0.70  Torr. 
Dots:  measured  distribution.  Curve:  calculated 
distribution. 
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VI.  DISCUSSION  AND  CONCLUSIONS 


The  Gaussian  velocity  distributions  used  to  calculate  the  results  shown 
in  Figs.  7-9  are  presented  in  Figs.  10,  11,  and  12,  together  with  the 
corresponding  Maxwellians  having  the  same  peak  velocity.  These  figures  reveal 
that  the  atomic  beam  velocity  distributions  are  indeed  nonthermal  and  much 
narrower  than  Maxwellians.  Because  the  distribution  at  0.7  Torr  is  broader 
than  the  distributions  at  0.2  and  0.9  Torr,  it  is  clear  that  we  are  not 
measuring  distributions  narrowed  by  gas-dynamics  effects.  Instead,  incomplete 
therraalization  of  the  hydrogen  atoms  must  be  the  cause  of  the  narrowing.  This 
interpretation  is  consistent  with  the  width  of  the  distribution  becoming  close 
to  the  Maxwellian  width  as  the  discharge  pressure  increases. 

The  peak  velocities  correspond  to  kinetic  energies  of  about  0.06  eV,  much 
smaller  than  the  approximately  2  eV  of  kinetic  energy  available  per  atom 
imnediately  after  dissociation.  Because  the  atoms  must  lose  the  excess  energy 
without  broadening  significantly  the '.r  energy  distribution,  the  main  energy- 
loss  mechanism  must  be  an  inelastic  process  with  a  relatively  large  loss  per 
collision.  Vibrational  excitation  of  the  background  hydrogen  molecules  (0.5^6 
eV  energy-loss/collision)  seems  to  be  a  good  candidate  for  such  a  process. 
Further  studies  in  this  area  are  required. 

It  is  well  known  that  the  velocity  passband  of  the  hexapole  or  quadrupole 
magnets  used  as  state  selectors  in  hydrogen  masers  is  rather  narrow  (about  251 
full  width  at  half  maximum).  Our  findings  indicate  that  it  should  be  possible 
to  design  the  state-selecting  magnet  having  a  velocity  passband  that  matches 
as  closely  as  possible  the  fairly  narrow  velocity  distribution  of  the  hydrogen 
atoms  effusing  out  of  the  dissociator.  Thus  the  efficiency  with  which  the 
maser  uses  its  hydrogen  supply  would  be  greatly  enhanced. 
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Fig.  10.  Model  Velocity  Distributions  Used  for  Fig.  7. 
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DISTRIBUTION  (10  s/m) 


Fig.  11.  Model  Velocity  Distributions  Used  for  Fig.  8. 
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Fig.  12.  Model  Velocity  Distributions  Used  for  Fig.  9. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security 
projects,  specializing  in  advanced  military  space  systems.  Providing  research  support,  the 
corporation’s  Laboratory  Operations  conducts  experimental  and  theoretical  investigations  that 
focus  on  the  application  of  scientific  and  technical  advances  to  such  systems.  Vital  to  the  success 
of  these  investigations  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay  current 
with  new  developments.  This  expertise  is  enhanced  by  a  research  program  aimed  at  dealing  with 
the  many  problems  associated  with  rapidly  evolving  space  systems.  Contributing  their  capabilities 
to  the  research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer 
and  flight  dynamics:  chemical  and  electric  propulsion,  propellant  chemistry,  chemical 
dynamics,  environmental  chemistry,  trace  detection:  spacecraft  structural  mechanics, 
contamination,  thermal  and  structural  control:  high  temperature  thermomechanics,  gas 
kinetics  and  radiation;  cw  and  pulsed  chemical  and  excimer  laser  development, 
including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control,  atmos¬ 
pheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmospheric 
optics.  light  scattering,  state-specific  chemical  reactions  and  radiative  signatures  of 
missile  plumes,  sensor  out-of-field-of-vicw  rejection,  applied  laser  spectroscopy,  laser 
chemistry,  laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency  stand¬ 
ards,  and  environmental  chemistry. 

Electronics  Research  laboratory:  Microelectronics,  solid-state  device  physics, 
compound  semiconductors,  radiation  hardening;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  microwave  semiconductor 
devices,  microwave/millimeter  wave  measurements,  diagnostics  and  radiometry,  micro¬ 
wave/millimeter  wave  thermionic  devices;  atomic  time  and  frequency  standards; 
antennas,  rf  systems,  electromagnetic  propagation  phenomena,  space  communication 
systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals,  alloys, 
ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  nondestructive 
evaluation,  component  failure  analysis  and  reliability;  fracture  mechanics  and  stress 
corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures 
as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnetosphcric,  auroral  and  cosmic  ray  physics, 
wave-particle  interactions,  magnetosphcric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using 
atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth's 
atmosphere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate 
radiations  on  space  systems;  space  instrumentation. 


